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We derive the low energy effective theory on a brane in six-dimensional chiral supergravity. The 
conical 3-brane singularities are resolved by introducing cylindrical codimension one 4-branes whose 
interiors are capped by a regular spacetime. The effective theory is described by the Brans-Dicke 
(BD) theory with the BD parameter given by wbd = 1/2. The BD field is originated from a modulus 
which is associated with the scaling symmetry of the system. If the dilaton potentials on the branes 
preserve the scaling symmetry, the scalar field has an exponential potential in the Einstein frame. 
We show that the time dependent solutions driven by the modulus in the four-dimensional effective 
theory can be lifted up to the six-dimensional exact solutions found in the literature. Based on the 
effective theory, we discuss a possible way to stabilize the modulus to recover standard cosmology 
1 and also study the implication for the cosmological constant problem. 
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I. INTRODUCTION 



Recently, much attention has been paid to six-dimensional supergravity P, 0, HJ . The most intriguing property 
of six-dimensional supergravity is that the four-dimensional spacetime is always Minskowski even in the presence of 
qj> branes with tension. A 3-brane with tension induces only a deficit angle in the six-dimensional spacetime and the 
tension does not curve the four-dimensional spacetime within the brane. This feature is called self-tuning and it may 
solve the cosmological constant problem [4j, |5|. This is the basis of the supersymmetric large extra-dimension (SLED) 
proposal [f|. 

There have been several objections to the idea of self-tuning 0, §]. The self-tuning relies on the classical scaling 
property of the model. The six-dimensional equations of motion are invariant under the constant rescaling gMN — ¥ 
C*~) i ^9mn and — > e^ _w , where gMN denotes the six-dimensional metric and cj> is the dilaton field. Then there is a 
■ modulus associated with this scaling property. Ref. [8 ; ] derived an effective potential for this modulus. This modulus is 
shown to have an exponential potential. Then there must be a fine-tuning of parameters to ensure that the potential 
vanishes in order to have a static solution. This is the reason why the static solution always has vanishing cosmological 
t-H ' constant. However, if this fine-tuning is broken, the modulus acquires a runaway potential and the four-dimensional 
spacetime becomes non-static. Non-static solutions in six-dimensional supergravity have been derived and they are 
supposed to correspond to the response of the bulk geometry to a change of tension of branes [§, Qil EH E2 . 
, However, it is difficult to deal with an arbitrary change of tension with a brane described by a pure conical 
singularity. This is because if we put matter on the brane other than cosmological constant, the metric diverges at the 
position of the brane. Rece ntly , it was suggested that we can regularize the brane by resolving it by a codimension 
one cylindrical 4-brane [H, Ell EE EH- This type of models may be regarded as a variation of Kaluza-Klein/hybrid 
brane world Once the brane becomes a codimension one object, it is possible to put arbitrary matter on the 

brane without having the divergence of the metric. Then it becomes possible to study the effect of the change of 
tension on the four-dimensional geometry on the brane. 

There is another interesting issue of whether it is possible to recover conventional cosmology at low energies in 
six-dimensional models. Recent works have shown that it is impossible to recover sensible cosmology if one derives 
cosmological solutions by considering a motion of branes in a given static bulk spacetime [l8|, Ell . It was concluded 
that the time-dependence of the bulk spacetime should be taken into account. 

In this paper, we derive a four-dimensional effective theory for the modulus in six-dimensional supergravity with 
resolved 4-branes by extending the analysis of Ref. [2(| which studied the low energy effective theory in the Einstein- 
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Maxwell theory [2l|. Arbitrary matter and potentials for the dilaton on 4-branes are allowed to exist. We use 
the gradient expansion technique to solve the six-dimensional geometry assuming that the deviation from the static 
solution is small [H, [23|]. The gradient expansion method has been applied to various types of braneworlds [24j |. 
Using this method, it is possible to solve the non-trivial dependence of the bulk geometry on the four-dimensional 
coordinates. By solving the effective four-dimensional equations, we can derive the time-dependent solutions and 
compare them with the exact six-dimensional time dependent solutions found in the literature [ToL [ill \l2^ . It is also 
possible to study whether we can reproduce sensible cosmology at low energies or not. We also study the possibility 
to stabilize the modulus using the potentials for the dilaton on the branes along the line of Ref . [25| . 

The paper is organized as follows. In section II, basic equations are summarized. In section III, we solve the 
six-dimensional equations of motion using the gradient expansion method. In section IV, the effective theory on the 
regularized branes is derived by imposing junction conditions. Then we derive time dependent cosmological solutions 
in the effective theory and compare them with the exact six-dimensional solutions. The possible way to stabilize the 
modulus is discussed. Section V is devoted to conclusions. 



II. BASIC EQUATIONS 



The relevant part of the supergravity action we consider is 

r4 11 #4 



5' 



2 2 v YJ 4 2L 2 



(1) 



where <f> is the dilaton, M is the fundamental scale of gravity, (d(j>) 2 :— g MN ' Bm^n^i F 2 := FmnF , and 
Fmn = 9m A N — 8nAm is the field strength of the gauge field A M . For the moment we are interested in solving the 
6D bulk equations of motion. In Sec. IIVI we will add two 4-branes (at positions y = y±) and Lj denotes the different 
bulk curvature scales on either sides of the branes, see Fig. [T] We start with the axisymmetric metric ansatz 

g M Ndx K Ux N =L]e 2X ^^+ee 2 ^y^- x ^f{y)d9 2 + (2) 

where capital Latin indices numerate the 6D coordinates while the Greek indices are restricted to the 4D coordinates. 
The evolution equations along the y-direction are given by 



1 



(3) 



where n v = e A v7 /Lj, K? is the extrinsic curvature of y — constant hypersurfaces, K is its 5D trace, 5 R? is the 5D 
Ricci tensor and 5 -D^ is the covariant derivative with respect to the 5D metric. Here, /t — [i and 9. The Hamiltonian 
constraint is 



1 



IR + K$K£ -K 2 = - w ( F v mF vM - -F 2 ) e-f - 2(n»^) a + W) 2 + 



1 



and the momentum constraints are 



b D> (Kf - 5£k) = Wi F !i MF vM n y e-t' + n*a 



where n„ 



The Maxwell equations are given by 

V M (e-*F MN ) = 0, 

where Vm is the covariant derivative with respect to the 6D metric. The dilaton equation of motion is 



A/, 



J—F 2 e -+- — = 
4Af 4 2L 2 



(4) 



(5) 



(6) 



(7) 
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III. GRADIENT EXPANSION APPROACH 

In this section we will use the gradient expansion method [22|, [H| to solve the 6D bulk equations. We assume that 
the length scale I is of the same order of Lj. The small expansion parameter is the ratio of the bulk curvature scale 
to the 4D intrinsic curvature scale, 

e = £ 2 \R\. 

We expand the various quantities as 

V = h fi „(x)+eh$(y,x)+--- , i> = V> (0) + + ■ ■ ■ , </> = 



(0) (i) 



(0) (1) 

KJ> =K h 6 +eKg 



(0) (1) 

F y d =F y e +eF y g 



As to the other quantities, we follow [20( and first assume 



b, = £ 1 /2 6 (l/2) 



1 /o <1/2) 



i In y 1 ' W 

=e 1 ' 2 F™ +■■■ , F liV =sF iiV 



(8) 



(9) 



(1/2) 



and then will show that all the 0(e 1 ^ 2 ) quantities in fact vanish. Since d^Ag ~ e 1 / 2 d y Ag, we have F^g = e 1 / 2 F^ e 
+ • • • . We will show that this 0(e 1 ^ 2 ) term in F^g also vanishes. The bulk energy-momentum tensor contains terms 
like F^\F V> ~ but these do not contribute to the low energy effective theory as they are higher order in the gradient 
expansion. The 5D Ricci tensor is given by 



5 R» = e- (R»[h] - V^i) - V^V v i> 



5p 6 



(10) 
(11) 



and 5 Rg = 0(e 3 ^ 2 ), where ip := tp^ — A. R" [h] and 2? M are respectively the Ricci tensor and the covariant derivative 
constructed from h^ v (x). 



A. Zeroth order equations 

The 9 component of the Maxwell equations at zeroth order reads 

d y (a l e-* i0)+ ^ a) F°y e ) =0, 
while the equation of motion for the dilaton at zeroth order is given by 

^fe^dy^) + \ $)V-* m -* W - \e^^ = 0. 

The (fiv) and (99) components of the evolution equations are given respectively by 

/ h(^) + ( 4 ^ + ¥ + ^ (0) ) ¥ 

1 ( 1 Y „2\-<f> W -2i>W _ 1 2A+0<°> 



4 \M 2 i 



v e e 



r 



4VM2£^V 4 



(12) 
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(15) 
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and the Hamiltonian constraint becomes 



4/ 



d v a \ ( dya f d y f 
f 



(jk&) e 2 ^ (0, ^ (0, +/(^) 2 - e 2 ^ (0, . 



The solutions for the above equations are obtained as 



1 



and 



V> (0) (2/,x) = $(2;)+CT(:r), ( %,z) = - lny - #(as), F y e= M 2 t\e^ )+i,i0) = M 2 t\e a{x \ 



(16) 



(17) 



(18) 



where /i and g are integration constants. The momentum constraint implies d^a = 0, and therefore a = constant. 

(1/2) 

This immediately leads to F fJ- g= and hence F^g = C(e 3 / 2 ). In the following, we put a = without loss of generality. 
The 6D metric at the zeroth order is given by 



g MN dx M dx N = e*^ 



L 2 ^ + e 2 fd9 2 



+ a 2 (y)h tll/ (x)dx^dx" . 



(19) 



Then we can see that <&(x) is associated with the scaling symmetry gMN — > e^gMN and — > e"^ - ". In fact, we will 
find that a solution for h^ u is given by = e*?^ if the brane preserves the scaling symmetry, where ?7 M „ denotes 
the 4D Minkowski metric. 



B. First order equations 



At first order, the (fiv) component of the evolution equations is given by 



V7 
L7 e 



-$/2 



(1) 



I + M 

y 2/ 



(i) (i) 



AL 2 



m 1 4 M ' 



(20) 



where 



J" 



1 / <°) 



M 4 



(i) 



FygFy + F ve Fy e )e 



yO 1 



1 (0) 



M 4 



(21) 



The 4D Ricci tensor R£ does not depend on y because it is computed from which is a function of x 11 only and 
the index is raised by /i p „. 

The 4D traceless part of Eq. (|20j) is found to be 



(22) 



where we defined 



J; -(1/4)^$ and 



(23) 



where 2? 2 <I> := h^V^V^ and (P$) 2 := h^V^V^. The general solution to the above equation is given by 



c;(i), 



(24) 



where the traceless tensor <C"(x) is the integration "constant" to be fixed by the boundary conditions. 
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The 4D trace part of the evolution equations is 



V7 - 



*/2 



y + 2f) K > + y Ke 



R-V 2 $- -(X>$) 2 



i — * 



(25) 



and the (00) component of the evolution equations is 



V7 - 



*/2 



y / 



2y 1 V y J 4 4Lj y V 



(26) 



The Hamiltonian constraint at first order reduces to 



R-V 2 <$>- -(£>$) 2 



Li 



% 2/ 



(i) o (i) 
^ +- itf 



1 e-* m , 2/e 



y ' £? » 



+ fl^« (27) 



The dilaton equation of motion at first order reads 



Now we define convenient quantities 



i 



/7„-*/2 / (1) (1) 

£- e —\K£ + Kj 
Li y \ 



J" := + J ^ 



(28) 



(29) 



and 



o (1) (1) If 



V> 1 



V 



M*PL Iy fJ 



Li y 



(30) 



The evolution equations for these variables can be derived using Eqs. (j25|) -(f28 |) . With some manipulation one arrives 
at 



d v (y 2 Vfj) = \e* ,2 yLi 
d y (y 2 VJlc) = \e* /2 yLi 



R + V 2 <P + -(£>$) 2 
i?-3X> 2 $- 7 -{V§) 2 



(31) 
(32) 



The two equations have the same structure as that of Eq. (|22[) . The general solution for each evolution equation 
contains one integration "constant" which will be determined by the boundary conditions. 
In terms of the above variables, the momentum constraint equations are simplified to 



V 



„ (e* /2 K;) - V, (e*' 2 K) + e"/ 2 JV^ = 0. 



(33) 



IV. JUNCTION CONDITIONS AND EFFECTIVE THEORY ON A REGULARIZED BRANE 

Our choice of parameters /i, q implies that /(y) vanishes at and ys- These points are conical singularities that 
are sourced by 3-branes. In order to accommodate usual matter on the branes we need to resolve these singularities. 
We will use the regularization scheme of [13l . |25| . The conical branes are replaced with cylindrical codimension-one 
branes at y = y± and their interiors are filled with regular caps. See figure [T] for a sketch of the model. The geometry 
of the caps and the central bulk is described by the 6D solutions found in the previous section, with different curvature 
scales L+ (L-) for the north (south) cap and Lq for the central bulk. 

The action of each brane is taken to be 



Si 



brane 



v(4) 



eA A )(a /i E - e^ A ) 



(34) 



FIG. 1: Schematic representation of the bulk spacetime with two regularized caps. 



where q^o is the induced metric on the 4-brane, V{4>) and U{<p) are the couplings to the dilaton, and £ m is the 
Lagrangian of usual matter localized on the brane. At this stage we assume that the brane matter £ m does not couple 
to the dilaton field. We introduce a Stueckelberg field E, which is obtained by integrating out the massive radial 
mode of a brane Higgs field. The equation of motion for E gives the gradient expansion form of the solution as [20] 



E(0, x)=n9 + c {0) {x) + ec {i) {x) + 

where n must be an integer because of the periodicity 6 ~ 6 + 2n. 
The jump conditions for the Maxwell field are 



[[n M F MN e-*]] = -eU{d N Y, - eA N ), 



while for the dilaton field we have 



(35) 



(36) 



(37) 



where [[F )] Vb := lim e ^ {F\y b +e ~ F\yb-<d- Here and hereafter in this section all the quantities are evaluated at the 
position of the brane under consideration. The Israel conditions are given by 



^ rji 

' A(tot) 



where 



^(tot) = -VSjl + U 



(%E - e^)(cTE - eA>) - -<^(%E - e^)(9 A E - eA x ) 



+ T- v 



(38) 



(39) 



and Ty represents the matter energy-momentum tensor. 



A. Zeroth order 



At zeroth order in the gradient expansion the junction conditions (|36|) - (|38[) are written as 



Maxwell: 
Dilaton: 
Israel [}iv) : 
Israel (00) : 



v7 ( °> <E>/2 

-i—F y e ye ' 
Li 

V_7 1 -*/2 
Li y 

VT ( 3 d y f 



3t/(°) ( 



n — eA 



(0) 



\2y 2/ 

v_7 £ e -*/2 
y 



l 

s -*/2 
1 



dV (o) i dU (0) e 



-<i> 



d0(°) 2 d^O) £ 2 / 
1 

: ~F 



(n-e^^ 
-<i> 



(0) 



.(0) 



(40) 
(41) 
(42) 
(43) 
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The above conditions relate several parameters with each other, and the detail of the parameter counting of the 
configuration is found in Rcf. 25]. In particular, the dilaton jump condition (|41[) and the Israel condition (|43[) imply 



y(0) dV (0) i e -*/[/(o) dU (0) 



(n — eA : 



(0) 



0. 



2 <ty(o) 2 jpf V 2 d0(o) 
The classical scaling symmetry is preserved by the special choice of the potentials @, [25j 



(44) 



(45) 



With these potentials the junction conditions (|4"0")) - (j4"3")) put no constraints on $(x) and Eq. (|4~4"|) is trivially satisfied. 
In this case the first order analysis will provide the equation of motion for &(x), as will be seen in the next subsection. 
In the following, we assume that at the zeroth order, the potentials are given by (|4"S"|) , that is, [/^°- ) (0' ') = u^e~^° ^ 2 
and (^(°) ) = e^ l0) ^ 2 . Then we expand the potentials as follows: 



v(<t>) = + + 

uq) = c/(°)(^)) + £ (c/W(^)) + ^w), 

where (cf)^) and U^ 1 '^ ') stand for the deviations from the zeroth order potentials. 



(46) 
(47) 



B. First order 

The 4D traceless part of the Israel conditions at first order is given by 

1 



[M] 



M A M> 



where := T" — (1/ '4)6 . The 4D trace part of the Israel conditions reduces to 
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(i) 


(i) " 
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AM 4 



T\ — 



M 4 



AV 



UW 1 

M 4 M 4 



dVW 1 dU^ e~* 



eA, 



where we defined 



and 



AV = V^H^) + lu^(^)^r(n-eA^) 



A:= C —(n-eA 



(or 



>) [eA^ + ^-e^)^ 1 ) 
Using the zeroth order junction conditions, Eq. (|49[) simply gives 

[[£]] 

The (00) component of the Israel conditions is 



4M 4 A M 4 



T 1 



A 



1 



M 4 M 4 M 4 

-Lav + uW^) 



^-'-^-(n-cAr 2 



d0(°) 2 d0(°) £ 2 / 



and the dilaton jump condition is 



M 4 



1 g/W 
' M 4 
1 d 



A 



M 4 £ 2 / 

d 2y(0) 1 d 2jj(0) g-* 



M 4 



d(/)( )2 2 #(0)2 £2J 

_rf> 

e 



(AF)--C/«(^°))^ r ^fn-eA( 0) ) 



n — eA 

2 



(0) 



(48) 



(49) 



(50) 



(51) 



(52) 



(53) 



(54) 
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Using the fact that the zeroth order potential have the scale invariant forms (|45|) . the above two conditions are 
combined to give 

m = ^Tf--L 4 ±(AV)- 1 ± i AV. (55) 
Therefore, the momentum constraints become 

Vv ( e */z T ; - e^AVS") = Qt/ - ±(AV) - ^AV^j e^V,<i>. (56) 

In terms of the energy-momentum tensor integrated along the ^-direction, 

T; := 2rfV7e*/ 2 T/, (57) 



this can be rewritten as 



T> v Tl - \T 6 e V^. (58) 



To fix the integration constants completely, we need the boundary conditions at the north and south poles. Near 
a pole with the coordinate y = y p , where p — {N, S 1 }, we have / ~ y — y p . In order for the evolution equations (I22p . 
(|3ip. and (|32f to be regular at the poles, we require 

K;, JC, J< \y-y P \ 1/2 ^0. (59) 

Now we can determine all the integration constants included in the general solutions for M.", K. and J. Since the 
structure of the evolution equations and boundary conditions are identical for these three variables, we summarize 
the procedure to fix the integration constants in Appendix [3 and here we focus on the resulting effective theory on 
the brane. 

Using Eqs. (j48|) and (|52|) together with the solution for and K. in terms of R and we end up with the effective 
equations 



(60) 



where the 4D gravitational couplings are defined as 



4 ■= 2 ^ 2 M 4 ' with l *= l \ Liydy, (61) 



vs 

2 



; denotes a covariant derivative with respect to the induced metric q^ v = a^_h p ^, R pv [q + ] is Ricci tensor computed 
from q+ v and the potential integrated along the ^-direction is defined as 

AV = 2n£y r fe' s ' /2 AV. (62) 
The first order equations for J give the equation of motion for 

(e% = f (t + x X - T? + 2± { AV + ) - 4AF+) + |^ (V - T/ + 2 A(aF") _ 4 AF") . (63) 
For simplicity let us ignore the matter energy-momentum tensor and the potential on the south brane: T ^ = 



AV = 0. In the absence of the (09) component of the energy momentum tensor on the north brane, the 4D effective 
equations can be deduced from the action 



S cS = J d A x^P 



-3- (%+] - c^bd*^*'") - AV + + Ct 



(64) 



with the Brans-Dicke parameter wbd = 1/2 (see also Appendix B of Ref. [HI]). 
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C. The exact time-dependent solutions in the 4D effective theory 



We now consider cosmological solutions in the 4D effective theory and compare them with the known solutions to 
the full 6D field equations [laCElCLl- 

Let us assume = AV = and Tf° = 0. We consider the case where the first order potential is scale invariant 
form. Then AV + oc e~*/ 2 and AV + — const. := A/re 2 , . We go to the Einstein frame defined by 

h iai = e s> q+ /> (65) 

and then the equations of motion become 

K\h]-\^R[h] = -Ae- 2 *^ + 2$|^l" - $| A $I^;, (66) 

= -Ae- 2 *, (67) 

where | stands for the covariant derivative with respect to the Einstein frame metric h^ v . 

Taking to be a flat Friedman-Robertson- Walker metric, h^dx^dx" = A 2 (t)(— dr 2 + dx 2 ), the equations of 
motion reduce to 

A" 9 1 / A' \ 2 1 1 

-y- = -AA 2 e~ 2 * - -& 2 , (±) = iAA 2 e- 2 * + -$' 2 , (68) 

and 

$" + 2 ^$' = AA 2 e- 2 *, (69) 

where ' := d/dr. For A = a solution of the above equations is 

A 2 {t) = A lT + A 2 , $(r) = ±V3hxA(r) + A 3 , (70) 

where Ai (i = 1, 2, 3) are integration constants. For A^Oa solution is 

A(r)=e CT , $(r) = lnA(r), (71) 

where C 2 = A/2. This indicates that g+ = e^rj^v which is expected from the scaling symmetry. 
The brane scale factor ab and the "radion" \& are given by 

a h = e-*' 2 A 1 tf = e*. (72) 

The solution ([70)) gives the same 4D observables (the brane scale factor and radion) as an exact 6D solution found by 
Copeland and Seto (equation (70) of [II]). In the same way the solution (f7Tj) reproduces the 4D quantities of their 
equation (76). (The latter solution was first found by Tolley et al. [I3|.) Thus we show that the solutions to the full 
6D equations are reproduced by our 4D effective theory on a regularized brane with the scale invariant potential and 
with or without additional "tension" A. 

At the zeroth order, the amplitudes i;' ' and of the potentials are fine-tuned. If we change the amplitude of 
the scale invariant potentials, which is equivalent to add a cosmological constant in the 4D effective theory, we get a 
runaway potential for $ and the 4D spacetime becomes non-static. 



D. Breaking the scale invariance 



Finally, let us consider the case where the first order potential breaks the scale invariance. As the BD parameter 
is given by 1/2, this model violates the constraints coming from the solar system experiments unless the BD scalar 
<I> is stabilized. It is suggested that the potential on a brane can naturally stabilize the modulus. For example, if we 
consider potentials Vi(<f>^) = v^e"^ and t r i(0 < -°- ) ) = u^e*^ , the effective potential is given by 

AT = 2^V7e*/ 2 (^ e ~ s * + WJyr eHt+1) ^ n - e ^ 0) ) 2 ) ■ ( ?3 ) 

As we saw in the previous subsection, if we take s = 1/2 and t — —1/2, AV is independent of <i>. However, in general, 
it is possible to have a potential with a minimum by choosing s,t,v^ and appropriately (25| . Then $ can be 
stabilized and general relativity (GR) is recovered. 
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V. 



CONCLUSIONS 



In this paper, we derived the low energy effective theory in the six-dimensional supergravity with resolved 4-branes. 
The gradient expansion method is used to solve the bulk geometry. The resultant effective theory is a Brans-Dicke 
theory with the Brans-Dicke parameter given by wbd = 1/2. If we choose the dilaton potentials on the branes so that 
they keep the scaling symmetry in the bulk and if we tune their amplitudes then there is no potential in the effective 
theory and the modulus is massless. Thus the static four-dimensional spacetime has vanishing cosmological constant. 
It is also possible to obtain time-dependent solutions due to the dynamics of the modulus field and we showed that 
they are identified with the six-dimensional exact time dependent solutions found in [l~2l ]. 

Even if the potentials preserve the scaling symmetry, it was found that there appears an effective cosmological 
constant in the four-dimensional effective theory by changing the amplitude of the potentials. Then in the Einstein 
frame, the modulus field acquires an exponential potential and the static solution is no longer allowed. Again, we 
showed that the cosmological solutions obtained in the effective theory can be identified with the six-dimensional 
exact time dependent solutions found in [llj, [HI, [12] • 

Our effective theory allows us to discuss cosmology with arbitrary matter on the brane. As the BD parameter 
is given by 1/2, it is impossible to reproduce realistic cosmology without stabilizing the modulus field. As it was 
suggested by Ref. (25[, it is easy to generate a potential for the modulus $ with a minimum by breaking the scaling 
symmetry from the dilaton potentials on the branes. Then it is possible to reproduce GR at low energies. However, 
once we stabilize the modulus, the cosmological constant on a brane curves the four-dimensional spacetime in the 
same way as in GR. 

Our result would indicate that it is possible to reproduce sensible cosmology in this six-dimensional supergravity 
model at low energies but it would be difficult to address the cosmological constant problem in this set-up. However, 
we should mention that our effective theory is valid only up to the energy scale determined by inverse of the size 
of extra-dimensions. This condition is roughly given by Hi* < 1 where H is the Hubble parameter. If we consider 
scales smaller than gravity becomes six-dimensional and from the table-top experiments, is smaller than a few 
fim. Then for H > 10~ 2 eV, our Universe becomes six-dimensional and it is impossible to use the four-dimensional 
effective theory. In order to address the behaviour of the universe at high energies, we should deal with time-dependent 
solutions directly in six-dimensional spacetime. This remains an open question. 

Finally, we briefly make a comment on the limit where the codimension one branes are shrunk to codimension two 
objects. Our effective theory shows no pathological behaviour in this limit as long as the four-dimensional energy- 
momentum tensor integrated along the ^-direction remains finite. However, in this limit, the first order extrinsic 
(i) 

curvature diverges and then the first order correction to the four-dimensional metric diverges. Then it is not 
clear whether there is a physical meaning in this limit. This is related to a deep issue of whether it is possible to put 
ordinary matter on codimension 2 objects [26] and we also leave this problem as an open question. 
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APPENDIX A: SOLVING THE BULK EVOLUTION EQUATIONS 



All of the key evolution equations in the main text have the form of 



d y y 2 y/fK(y,x) 



e^/iyLjRix), 



(Al) 



subject to the boundary conditions 



K(y N ,x) = K(y Sl x) = 0, 



(A2) 
(A3) 
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In the south cap we have the solution 



K= vl_yjL e */2 LR 



In the bulk the solution can be written as 



K 



e*/ 2 [y 2 L R + C(x)] 



where the integration constant C(x) is determined by the condition (|A3p as 

C{x) = [-y 2 _L + (y 2 _ y%)L_] R + 2y 2 v T< "■' " ' . 
In the north cap we have the solution 

K = ll-Ae^L+R, 

and the boundary condition (|A3|) requires 

[(y 2 N - y 2 + )L+ + y\ L ] R + C = -2y\^Ue^' 2 T + . 
In the above we defined f± := f(y±). Substituting Eq. (|A6|) into Eq. (|A8|) we obtain 



rVN \ 

/ L 7 ydy ■ i? = - Vi V^e-* /2 T 
■>j/s / i= ± 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



APPENDIX B: 0{e L/I ) QUANTITIES 



The /i component of the 0(e 1 ^ 2 ) Maxwell equations reads 

(1/2) \ 



and thus we have 



d y [ e~+ K V JW»J = 0, 



F"= M 2 - 



(Bl) 



(B2) 



The ©(e 1 / 2 ) evolution equation reduces to 

e -«/2 i 



I (0) 11/2) 

M 4 V 
2T 



which can be integrated to give 



(1/2) e 



*/2 



y 2 V7 



v 



(1/2) (1/2) 



(B3) 



(B4) 



The 0(e) evolution equations contain terms like F^yF^tx h^\C^C\ / f(y)- In the cap regions, we thus require 



(1/2) 



C( = to avoid the singular behavior at the poles. Further, the regularity of Kg at the poles imposes = in 
the cap regions. To fix the integration constants in the central bulk, we use the Maxwell jump conditions and Israel 
conditions at each branc: 



n y F vfi e 



(1/2) 



-eU (d^-cA^) 



U 



(1/2) 



(B5) 
(B6) 
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(1/2) (1/2) 

Combining these two equations and noting that F ytl = =Kg in each cap, we obtain two linear algebraic equations 
for the bulk values of and : 



(1/2) | 

K," 



eM 4 



(n - cAf^ ,,,, 



(1/2) 



= 0. 



Therefore, C( = C% = in the bulk. Now we also see that 

- eA^) (1/2) = on the branes, 

and 6 M =C(e 3 / 2 ). 



(B7) 



(B8) 
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